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SUMMARY 


An  exact  theoretical  method  is  developed  that  permits  the  determina¬ 
tion  of  the  pressure  field  of  a  wing-body  combination  having  a  circular 
body  and  a  wing  with  supersonic  leading  and  trailing  edges •  Detailed 
calculations  have  been  performed  for  wing-body  combinations  composed  of 
rectangular  wings  mounted  at  incidence  on  bodies  at  zero  angle  of  attack 
for  effective  chord-radius  ratios  of  4  or  less.  For  large  effective 
chord-radius  ratios  some  asymptotic  results  have  been  obtained*  It  was 
determined  that  for  the  family  of  combinations  having  an  effective  chord- 
radius  ratio  of  4  the  area  of  the  wing  blanketed  by  the  body  does  not 
generate  any  lift  itself  but  rather  acts  to  support  the  lift  generated 
by  the  exposed  wing,  and  that  the  body  is  less  than  50  percent  effective 
in  reflecting  lift  back  to  the  wing.  For  chord-radius  ratios  less  than  4 
the  relative  amount  of  reflection  increases. 

The  significant  fact  was  determined  that  for  rectangular  wing-body 
combinations  for  which  the  effective  chord-radius  ratio  is  greater  than  4 
most  of  the  loss  of  lift  due  to  interference  can  be  estimated  from  the 
first  term  of  the  Fourier  series  used  in  the  analysis.  This  fact  was 
used  to  determine  asymptotic  lift  results  for  the  region  where  no  exact 
calculations  were  made.  The  asymptotic  expressions,  together  with  the 
calculations,  allowed  the  construction  of  design  charts  showing  the  lift 
and  c enter -of-pres sure  location  of  the  exposed  wing  panels  as  a  function 
of  effective  aspect  ratio  and  effective  chord-radius  ratio.  The  charts 
show  that  as  a  result  of  interference  the  lift  on  the  exposed  wing  panels 
in  combination  with  the  body  can  be  reduced  as  much  as  15  percent  below 
the  value  for  the  wing  panels  joined  together,  and  that  the  center  of 
pressure  of  the  exposed  wing  panels  can  move  forward  as  much  as  4  percent 
of  the  wing  chord. 


INTRODUCTION 


In  recent  years  the  problems  of  supersonic  wing-body  interference 
have  occupied  the  attention  of  many  workers  in  aerodynamics.  The  large 
amount  of  effort  expended  on  the  subject  is  a  result  of  the  important 
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effects  that  interference  can  have  on  the  over-all  aerodynamic  character¬ 
istics  of  wing-body  combinations.  The  trend  toward  using  large  bodies 
and  small  wings  at  supersonic  speeds,  especially  for  missiles,  is  the 
prime  reason  for  the  increased  importance  of  wing-body  interference  at 
these  speeds. 

Much  significant  work  has  already  been  done  in  the  field.  In  ref¬ 
erence  1,  Spreiter  has  shown  that,  when  a  wing-body  combination  is  slender 
in  the  sense  of  his  paper,  simple  expressions  for  the  lift  and  moment 
coefficients  can  be  derived.  These  results  were  obtained  by  reducing  a 
three-dimensional  problem  for  the  wave  equation  to  a  two-dimensional 
problem  for  Laplace Ts  equation.  Another  approach  is  that  of  simplifying 
the  differential  equation  by  using  conical  boundaries.  Following  this 
approach,  Browne,  Friedman,  and  Hodes  in  reference  2  obtained  a  solution 
for  the  pressure  field  of  a  wing-body  combination  composed  of  a  flat 
triangular  wing  and  a  cone  both  with  a  common  apex.  The  use  of  all- 
conical  boundaries  reduces  the  problem  to  one  of  conical  flow  for  which 
powerful  methods  of  solution  are  available. 

Several  investigators  have  presented  methods  for  determining  the 
pressure  field,  including  the  effect  of  interference,  acting  on  wing-body 
combinations  employing  circular  fuselages  and  wings  not  necessarily 
slender.  In  reference  3>  Ferrari  has  given  an  approximate  method  of 
obtaining  the  "interference  of  the  wing  on  the  streamlined  body,  assuming 
that  the  induced  field  generated  by  the  wing  is  that  which  would  exist 
around  the  wing  if  it  were  placed  in  the  uniform  stream  alone."  Simi¬ 
larly,  the  interference  of  the  body  on  the  wing  has  been  determined.  The 
results  of  Ferrari  thus  represent  a  first  approximation  and,  while  a 
second  approximation  using  the  method  is  possible  in  principle,  it  appears 
that  too  much  labor  would  be  involved. 

Another  method  for  estimating  the  effect  of  interference  on  the  aero¬ 
dynamic  properties  of  wing-body  combinations  which  are  not  necessarily 
slender  is  given  in  reference  k.  In  this  reference  the  method  is  applied 
to  determining  the  drag  of  symmetrical  wing -body  combinations;  it  is  also 
applicable  to  the  calculations  of  the  lifting  pressures  acting  on  combi¬ 
nations  employing  wings  with  supersonic  edges.  In  reference  an  essen¬ 
tially  new  method  of  solving  a  wide  class  of  wing-body  interference  pro¬ 
blems  has  been  presented.  The  method  is  based  on  decomposing  the  inter¬ 
ference  of  a  wing-body  combination  into  a  number  of  Fourier  components 
and  solving  the  problem  for  each  component  in  a  manner  similar  to  that 
used  by  von  Karman  and  Moore  in  reference  6  for  bodies  of  revolution. 

To  summarize  the  present  situation,  it  can  be  said  that  the  existing 
interference  methods  fall  into  two  categories.  The  first  category 
includes  those  attempts  to  solve  difficult  boundary -value  problems  as  in 
references  3>  5,  and  7*  These  methods  have  the  shortcomings  that  they 

all  involve  a  great  deal  of  labor,  and  that,  with  the  exception  of  ref¬ 
erence  5,  they  represent  approximations  to  the  true  linear-theory  solu¬ 
tions.  The  second  category  includes  simpler  methods  such  as  those  of 
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references  1,  2 ,  and  8  and  existing  approximate  engineering  methods. 

TOiese  methods  have  the  disadvantages  of  being  applicable  only  to  special¬ 
ized  configurations  or  of  being  based  on  assumptions  of  unknown  validity. 

At  the  present  time,  there  is  a  definite  need  for  a  simplified  gen¬ 
eral  theory  of  wing-body  interference  that  will  give  results  of  engineer¬ 
ing  accuracy  for  a  wide  range  of  wing  aspect  ratio,  taper  ratio,  leading - 
edge  sweepback  angle,  and  span -diameter  ratio.  The  reasons  that  none  has 
yet  been  developed  are  twofold:  first,  the  lack  of  systematic  experi¬ 
mental  results  makes  an  empirical  theory  difficulty  and  second,  incomplete 
knowledge  of  the  mechanism  of  wing-body  interference,  because  of  the 
shortage  of  exact  solutions,  makes  it  difficult  to  develop  a  reliable 
rational  method.  The  development  of  a  simple  general  theory  of  wing-body 
interference  seems  thus  to  be  dependent  on  obtaining  a  few  exact  solutions 
that  will  give  sufficient  insight  into  the  mechanism  of  interference  to 
permit  valid  simplifying  assumptions.  Such  solutions  are  also  necessary 
to  assess  the  validity  of  the  assumptions  underlying  present  theories. 

In  reference  5>  3-  solution  has  been  obtained  to  a  wing -body  problem 
that  is  exact  within  the  limitations  of .linear  theory.  It  is  the  pur¬ 
pose  of  this  paper  to  present  a  resume  of  reference  5  with  special 
emphasis  on  the  physical  principles  underlying  the  interference.  It  is 
hoped  that  the  insight  gained  into  the  mechanism  of  interference  in  this 

paper  will  be  sufficient  to  allow  the  development  of  a  simple,  reliable 
engineering  method. 

'khis  paper,  a  physical  description  of  the  interference  problem 
w  11  be  given,  and  an  illustrative  example  of  a  rectangular  wing  at  inci¬ 
dence  on  a  round  body  at  zero  angle  of  attack  will  be  presented  in  detail. 


SYMBOLS 


a  body  radius 

A  aspect  ratio  of  wing  formed  by  joining  exposed  half -wings 

together 

c  chord  of  rectangular  wing 

c*  effective  chord-radius  ratio 

CD  minimum  wave  drag  coefficient  of  wing  based  on  exposed 

wing  area 

cr  chord  at  wing-body  juncture 

ct  chord  at  wing  tip 
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f2nW 

g2n(x) 

h2n(x) 


"W 


^2n(T>r ) 


JW 


hlC 


2n 


M2n(x )  ^2n(x ) 


n 


r2Xi 

q. 

r, 0,x 


velocity  amplitude  function  of  n!th  Fourier  component 

axial  strength  function  of  nTth  Fourier  component 

axial  strength  function  for  unit  step  in  f2n(x) 
origin 

Lwc 

% 

function  used  in  determining  P2n 

lift  on  exposed  half-wings  joined  together 

lift  on  exposed  half -wings  in  combination  with  body 

lift 

lift  due  to  n!th  Fourier  component 
odd  integer 

free-stream  Mach  number 
characteristic  functions 
number  of  Fourier  component 

pressure  coefficient 

interference  pressure  coefficient  due  to  nTth  Fourier 
component 

free-stream  dynamic  pressure 

cylindrical  coordinates:  y  =  r  cos  0,  z  =  r  sin  6 
(See  fig.  1. ) 

semispan  of  wing-body  combination 
maximum  wing  thickness 

thickness  ratio  of  double-wedge  wing 


u,v,w 


axial,  lateral,  and  vertical  perturbation  velocities, 
respectively 


MCA  TN  2677 


x,y,z 


xc  .p. 


V 

Og 

aZ 

a-u 

3 

3A 

& 

X 


t>T 


A 

9 

9C 

'Pan 

92n*(r,e,x) 

t(x,r) 


Cartesian  coordinates:  x,  axial  coordinate;  y,  lateral 
'  coordinate;  z,  vertical  coordinate 
(See  fig.  1.) 

center  of  pressure  of  wing  lift  measured  from  wing  lead 
ing  edge 

free-stream  velocity 
body  angle  of  attack 
local  angle  of  attack 

upwash  angle  of  body-alone  flow 
wing  angle  of  attack 

Vm2-i 

effective  aspect  ratio 
sin”1x 
taper 

dummy  variables  of  integration 
sweep  angle  of  wing  leading  edge 
interference  perturbation  velocity  potential 
combination  perturbation  velocity  potential 
wing-alone  perturbation  velocity  potential 
n’th  Fourier  component  perturbation  velocity  potential 
fundamental  solution  for  n’th  Fourier  component 
undetermined  function  of  x  and  r 

PHYSICAL  PRINCIPLES 


Prior  to  a  mathematical  formulation  of  the  wing-body  interference 
problem,  it  is  well  to  define  interference  and  to  explain  how  it  arises 
With  a  stationary  wing  or  a  stationary  body  in  a  uniform  parallel  flow. 
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there  are  associated  the  wing-alone  and  body-alone  flow  fields*  If  the 
wing  is  immersed  in  the  body-alone  flow  field,  the  fluid  velocity  due  to 
the  body-alone  field  will  not,  in  general,  be  tangential  to  the  wing  (or 
conversely).  For  this  reason,  the  sum  of  the  body -alone  plus  wing-alone 
flow  fields  will  not  be  the  flow  field  for  the  body  and  wing  together. 

The  difference  between  the  flow  field  of  the  body  and  wing  together  and 
the  sum  of  the  body-alone  and  wing-alone  flow  fields  is  defined  to  be  the 
interference  flow  field. 


Effect  of  Forebody  on  Wing 


The  effects  of  wing-body  interference  on  the  flow  field  of  a  wing- 
body  combination  are  illustrated  by  considering  separately  the  effects  of 
each  component  on  the  others.  For  the  purposes  of  this  discussion 
figure  1  shows  a  wing-body  combination  divided  into  the  part  in  front  of 
the  leading  edge  of  the  wing-body  juncture,  henceforth  called  the  nose, 
the  winged  part  and  the  part  behind  the  wing  trailing  edge,  henceforth 
called  the  afterbody.  The  wing  may  be  twisted  and  cambered,  but  it  is 
assumed  to  have  supersonic  leading  edges.  Consider  now  the  flow  as  it 
progresses  past  the  body.  At  the  body  nose  the  flow  is  that  around  a 
body  of  revolution,  and  it  can  be  treated  by  existing  methods  such  as 
those  of  references  1,  6,  and  9*  When  the  body  is  at  angle  of  attack  0,3, 
there  is  an  upwash  field  in  the  horizontal  plane  of  symmetry  of  the  body. 
If  the  body  is  sufficiently  slender,  the  flow  field  in  a  plane  at  right 
angles  to  the  body  axis  corresponds  to  that  around  a  circular  cylinder 
In  a  uniform  stream  of  velocity,  V  sin  0,3.  This  gives  an  upwash  field 
in  the  horizontal  plane  of  symmetry  of  the  body  of 

Ou  =  <xB  (1  +  a2/y2)  (1) 

The  effect  of  this  upwash  on  the  wing  can  be  obtained  by  considering  the 
wing  to  be  at  angle  of  attack  and  twisted  according  to  equation  (l)  and 
by  applying  the  formulas  of  supersonic  wing  theory.  The  wing  pressure 
field  so  obtained  is  exact,  within  the  limitations  of  the  theory,  for 
that  section  of  the  wing  outboard  of  the  Mach  line  emanating  from  the 
leading  edge  of  the  wing-body  juncture.  If  the  wing  is  located  close  to 
the  body  nose  so  that  there  is  a  chordwise  variation  in  the  upwash  field 
due  to  the  body,  then  the  wing  is  effectively  cambered,  and  the  solution 
is  more  difficult.  However,  for  most  wing-body  combinations  It  is  pos¬ 
sible  to  disregard  the  effect  of  the  nose,  and  to  assume  that  the  wing 
is  attached  to  a  circular  cylinder  that  extends  upstream  indefinitely. 


Mutual  Effects  Between  Body  and  Wing 


The  mutual  interference  between  the  body  and  wing  on  the  winged  part 
of  a  combination  causes  an  interference  field  acting  on  the  body  and  on 
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the  wing  inboard  of  the  Mach  line  emanating  from  the  leading  edge  of  the 
wing -body  juncture. 

The  wing-alone  flow  field  does  not,  in  general,  produce  flow  tan¬ 
gential  to  the  position  to  be  occupied  by  the  body  surface.  An  inter¬ 
ference  flow  field  must  arise  that  cancels  the  velocity  induced  by  the 
wing-alone  flow  field  normal  to  the  body  while  not  changing  the  body 
shape.  Alternately,  the  origin  of  the  interference  field  can  be  explained 
in  the  following  manner  •  The  wing  and  body  can  be  thought  of  as  sources 
of  pressure  disturbances  that  radiate  in  all  directions  in  downstream 
Mach  cones.  The  wing  disturbances  which  radiate  toward  the  body  are  in 
part  reflected  back  by  the  body  onto  the  wing  and  in  part  transmitted 
onto  the  body  giving  rise  to  interference  pressures.  Likewise,  the  dis¬ 
turbances  originating  on  the  body  pass  onto  the  wing  and  affect  the  pres¬ 
sures  there.  It  is  apparent  that  the  determination  of  the  interference 
pressure  field  on  the  body  and  on  the  wing  inboard  of  the  Mach  line  of 
the  juncture  is  the  crux  of  the  wing -body  interference  problem. 


Mutual  Effects  Between  Wing  Panels 


To  determine  the  region  of  influence  of  one  wing  panel  on  another, 
it  is  necessary  to  trace  the  path  of  a  pulse  from  one  wing  panel  across 
the  body  onto  the  other.  The  path  traced  across  the  body  by  the  pulse 
originating  at  the  leading  edge  of  the  wing-body  juncture  is  the  forward 
boundary  of  the  region  of  influence  of  one  wing  panel  on  the  body.  (See 
fig.  1.)  It  is  clearly  the  helix  intersecting  all  parallel  elements  of 
the  cylinder  at^the  Mach  angle.  The  boundary  crosses  the  top  of  the  body 

a  distance  of  — ^/m2-1  downstream  and  reaches  the  opposite  wing-body 

juncture  a  distance  ira  M2-l  downstream.  A  pulse  originating  at  a 
,  point  on  one  wing  panel  and  traveling  to  a  point  on  the  other  panel  can 
travel  around  the  body  on  its  surface  to  the  opposite  juncture  and  then 
along  the  wing  to  a  given  point,  or  it  can  leave  the  body  tangentially 
before  reaching  the  opposite  wing  juncture  in  a  straight  path  to  the 
point.  The  second  means  of  transmitting  the  impulse  Is  shorter  in  dis¬ 
tance  than  the  first  and  is  the  one  which  determines  the  forward  bound¬ 
ary  of  the  region  of  influence  of  one  wing  panel  on  the  other.  Applying 
this  consideration  to  the  pulse  originating  at  the  leading  edge  of  one 
wing -body  juncture,  it  is  easy  to  show  that  the  forward  boundary  of  the 
region  of  influence  of  one  wing  panel  on  the  opposite  wing  panel  is 
given  by  the  equation 


This  boundary  is  also  shown  in  figure  1,  and  it  becomes  parallel  to  the 
Mach  direction  at  distances  far  from  the  body. 
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Effects  on  the  Afterbody 

fir 


As  far  as  the  interference  effect  of  the  body  on  the  wing  is  con¬ 
cerned,  it  is  confined  to  the  winged  part  of  the  combination,  but  the 
effect  of  the  wing  on  the  body  is  felt  also  on  the  afterbody.  If  the 
downwash  were  known  everywhere  in  the  wing  wake,  then  the  wake  could  be 
considered  as  an  extension  of  the  wing  with  twist  and  camber.  The  wing 
wake  and  afterbody  could  then  be  incorporated  with  the  winged  part  of 
the  combination  and.  treated  in  the  same  manner.  However,  the  actual 
downwash  pattern  in  the  wing  wake  depends  on  the  interference  effect  of 
the  body  on  the  wing.  It  is  thus  apparent  that  the  solution  of  the  after¬ 
body  problem  requires  that  the  interference  problem  for  the  winged  part 
of  the  combination  be  solved  first.  Only  the  winged  part  of  the  combina¬ 
tion  is  analyzed  in  detail  in  this  report. 


MATHEMATICAL  FORMULATION  OF  THE  PROBLEM 


Throughout  the  analysis,  the  body  radius  is  taken  as  unity  and  M2 
is  taken  as  2  so  that  J3  =  1.  Any  formula  can  be  generalized  to  any  body 
radius  by  dividing  all  length  symbols  by  a,  and  to  any  Mach  number  by 
dividing  all  streamwise  lengths  by  p,  by  multiplying  all  pressure  and 
lift  coefficients  by  p,  and  leaving  all  potentials,  lift  forces,  and 
span  loading  unaltered.  It  is  necessary  to  specify  the  wing  alone  before 
any  detailed  interference  calculation  can  be  carried  out.  However,  in 
the  theoretical  solution  of  the  problem  the  wing-alone  definition  is 
arbitrary.  The  flow  field  about  the  combination  does  not  depend  on  the 
definition  of  the  wing  alone. 

The  analysis  is  confined  to  the  cases  for  which  the  aspect  ratios 
are  sufficiently  large  that  the  Mach  lines  from  the  leading  edges  of  the 
wing  tips  do  not  intersect  the  wing -body  juncture.  This  is  the  case  if 
the  aspect  ratio  obeys  the  following  inequality: 


Under  such  circumstances,  the  tips  have  no  influence  on  the  wing -body 
interference  of  the  winged  part  of  the  combination. 

The  wing-body  combinations  being  considered  in  this  report  are  those 
possessing  a  vertical  plane  of  symmetry.  To  eliminate  any  questions  of 
the  symmetry  of  the  boundary  conditions  above  and  below  the  wing,  the 
assumption  is  made  that  the  wings  have  supersonic  leading  and  trailing 
edges.  This  assumption,  together  with  that  concerning  tip  effects,  is 
sufficient  to  insure  the  fact  that  pressure  communication  between  the 
upper  and  lower  wing  surfaces  can  have  no  effect  on  the  interference  on 
the  winged  part  of  the  combination.  Thus  it  is  immaterial  whether  the 
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problem  at  hand  is  symmetrical  or  antisymmetrical  with  regard  to  the 
horizontal  plane  of  symmetry.  Attention  will  henceforth  be  focused  above 
the  horizontal  plane  of  symmetry. 

The  assumptions  with  regard  to  tip  effects  and  the  sweep  of  leading 
and  trailing  edges  enable  the  drag  of  a  symmetrical  wing-body  combination 
at  zero  angle  of  attack  to  be  obtained  from  solutions  for  lifting  sur¬ 
faces.  Suppose  the  problem  solved  is  that  for  a  flat  wing  at  angle  of 
attack  aw  mounted  on  a  body  at  zero  angle  of  attack.  The  pressures 
are  antisymmetric  with  respect  to  the  horizontal  plane  of  symmetry.  How¬ 
ever,  considering  the  same  pressures  as  positive  and  symmetrical  yields 
the  solution  for  the  pressure  field  of  a  single-wedge  wing  of  half 
angle  aw  at  zero  angle  of  attack.  This  fact  is  of  use  in  obtaining  the 
drag  coefficient  at  zero  angle  of  attack,  as  will  be  illustrated. 


General  Decomposition  of  Boundary-Value  Problem 


The  general  case  of  the  wing  and  body  at  different  angles  of  attack 
is  considered.  Following  the  suggestion  of  Lagerstrom  and  Van  Dyke  in 
reference  10,  the  problem  can  be  broken  down  into  several  simpler  prob¬ 
lems  on  the  assumption  that  the  wing  boundary  conditions  can  be  applied 
in  the  plane  of  the  wing  and  those  for  the  body  can  be  applied  on  a 
cylindrical  surface.  Figure  2  shows  the  decomposition  of  the  wing-body 
combination  into  three  components.  Component  (a)  is  simply  the  body 
alone,  which  creates  an  upwash  field  in  the  region  that  is  to  be 

occupied  by  the  wing.  Components  (b)  and  (c)  are  combinations  with  wings 
of  the  same  plan  form;  but  while  component  (b)  has  the  wing-alone  angle 
of  attack,  component  (c)  has  a  wing  with  angle  of  attack  -a^.  The  sig- 
n^:*‘Gance  "this  particular  method  of  decomposing  the  general  wing-body 
problem  is  that  component  (a),  the  body  alone,  can  be  solved  by  known 
methods  and  components  (b )  and  (c)  with  bodies  at  zero  angle  of  attack 
can  be  solved  by  the  methods  of  this  report.  Henceforth  the  analysis  is 
confined  to  combinations  with  bodies  at  zero  angle  of  attack. 


Wing  alone  and  interference  potentials.-  Consider  now  a  combination 
with  the  body  at  zero  angle  of  attack  and  let  cpp  be  its  potential. 

(See  fig.  3.)  This  potential  can  be  considered  the  sum  of  a  wing-alone 
potential  cp^-  and  of  an  interference  potential  cp. 


<P„  = 


% 


+  <p 


u 


The  essential  problem  is  to  determine  cp-  First,  select  a  convenient 
way  of  extending  the  wing  through  the  body  to  form  the  wing  alone,  thereby 
specifying  Ty.  The  wing-alone  flow  field  in  general  produces  veloc- 

ities  normal  to  the  surface  that  will  enclose  the  circular  cylinder 

as  illustrated  in  figure  3  for  the  region  above  the  wing.  In  figure  3 
and  the  subsequent  figure,  all  bodies  are  shown  as  cylinders  parallel 
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to  the  x  axis.  While  the  bodies  of  the  actual  configurations  in  some 
cases  are  slightly  distorted  cylinders,  they  are  nevertheless  shown  as 
true  cylinders.  This  procedure  is  compatible  with  the  fact  that  the 
boundary  conditions  are  to  be  applied  on  a  true  cylinder.  The  value 

of  varies  with  6  and  with  x.  This  means  that  a  body  conforming 

dr 

to  the  wing-alone  flow  field  is  distorted  in  a  complicated  fashion.  Wow 
since  the  body  must  be  circular,  there  must  arise  an  interference  poten¬ 
tial  cp  that  identically  cancels  at  the  body  surface,  thereby 

or 


straightening  it. 


at  r  =  1  (4) 

or  or 


There  are  two  other  conditions  to  be  fulfilled  by  9*  It  must  not  dis¬ 
tort  the  shape  of  the  wing  when  added  to  9y  to  produce  9q-  Thus 
when  0  =  0 , 


(5) 


or  cty  =  0  for  the  interference  combinations  as  shown  in  figure  3-  The 
last  condition  is  that  the  interference  potential  must  be  zero  ahead  of 
the  winged  part  of  the  combination. 

cp=0  x<0  (6) 


Equations  (4),  (5)>  and  (6)  are  the  essential  boundary  conditions  on  9. 

The  normal  velocity  —  to  be  induced  at  the  body  surface  by  the 

dr 

interference  potential  can  be  analyzed  at  any  given  streamwise  position 
as  a  Fourier  cosine  series.  The  amplitudes  of  the  various  Fourier 
cosine  terms,  f2n(x),  vary  with  x,  the  streamwise  distance.  Thus, 


•* 


if 


~  ^  f2n(x)  cos  2n@  =  “  at  r  =  1  (7) 

n=o 

Only  even  multiples  of  0  are  considered  because  of  the  vertical  plane 
of  symmetry.  Consider  that  the  interference  potential  is  decomposed  into 
a  series  of  potentials  such  that  each  cancels  one  Fourier  component  of 
the  velocity  at  the  body  surface,  that  is, 

00 

*  ■  I 


n=o 


(8) 


NACA  TN  2677 


11 


with 


frPan 

Sr 


=  f2n(x)  cos  2n0 


at 


r  =  1 


(9) 


Then  the  combination  giving  the  interference  potential  cp  can 
posed,  in  a  series  of  combinations  each  giving  one  of  the  cp^ 
The  decomposition  is  illustrated  in  figure  4.  For  n  =  0, 


be  de corn- 
values  . 


and  there  is  no  variation  of  the  normal  velocity,  pressure,  or  potential 
with  0.  Thus  the  first  interference  combination  is  a  body  of  revolution. 
The  pressure  field  acting  on  the  body  of  such  a  combination  can  be  deter¬ 
mined  by  the  method  of  reference  11.  This  n  =  0  interference  combina¬ 
tion  has  the  very  simple  significance  that  its  flow  normal  to  the 

r  =  1  cylinder,  ~°,  subtracted  from  ^2  reduces  the  flow  across  the 

or  Sr 

body  to  zero  when  averaged  from  0  =  0  to  0  =  or  at  any  streamwise  loca¬ 
tion.  For  n  =  1, 


Sr 


=  f2(x)  cos  20 


and  the  normal  velocity,  pressure,  and  potential  will  vary  as  cos  20. 

To  summarize  briefly,  it  has  been  shown  that  the  general  inter¬ 
ference  problem  of  a  body  and  wing  at  different  angles  of  attack  can  be 
broken  down  into  wing-body  problems  with  bodies  at  zero  angle  of  attack 
as  shown  in  figure  2.  Combinations  with  the  body  at  zero  angle  of  attack 
are  decomposed  into  wings  alone  plus  interference  combinations  as  in 
figure  3-  The  interference  combinations  are  finally  decomposed  into 
their  Fourier  components  as  in  figure  4. 

A  general  method  for  determining  the  characteristics  of  any  Fourier 
component  will  now  be  given.  It  will  be  shown  that  good  accuracy  can  be 
obtained  for  the  interference  potential  with  few  Fourier  components. 


Method  of  Solution 


The  problem  to  be  solved  is  that  of  a  supersonic  wing-body  combina¬ 
tion  subject  to  the  conditions  already  mentioned,  but  with  the  wing  and 
body  possibly  at  different  angles  of  attack.  This  problem  is  reduced  to 
a  body-alone  problem  and  two  wing-body  problems  with  the  body  at  zero 
angle  of  attack  as  shown  in  figure  2.  The  body-alone  problem  can  be 
solved  by  existing  methods  such  as  references  6  and  The  procedure 
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necessary  to  solve  either  wing-body  problem  as  given  in  reference  5  is 
now  summarized. 


Wing-alone  potential  and  velocity  amplitudes.-  The  first  step  in 
the  solution  is  to  determine  the  wing-alone  potential,  cp^.  The  methods 
of  wing  theory  available  for  doing  this  will  not  be  discussed  here. 

From  the  value  of  Ty  the  value  of  the  normal  velocity  induced  at 
the  r  =  1  cylinder  follows  readily.  This  normal  velocity  distribution 
is  expanded  in  a  Fourier  cosine  series  of  even  multiples  of  6 . 


Sr 


00 

f2n(x)  cos  2n 6  at  r  =  1 

n=o 


(10) 


The  functions  fPn(x)  are  called  the  velocity  amplitude  functions  and 
are  part  of  the  boundary  conditions  for  the  interference  potential. 


Interference  potential.-  The  solution  for  the  interference  potential 
is  obtained  by  summing  the  potentials  for  the  series  of  Fourier  component 
potentials.  The  determination  of  the  potential  for  any  one  Fourier  com¬ 
ponent  (as  illustrated  by  fig.  k)  has  been  carried  out  formally  in  ref¬ 
erence  5  by  Laplace  transform  theory.  It  was  shown  that  the  formal 
mathematical  solution  is  equivalent  to  distributing  fundamental  solutions 
for  cp  along  the  axis  of  the  combination  in  a  manner  similar  to  that  of 
von  Karman  and  Moore  in  reference  6.  Let  cp2n*(x-£,r,0)  be  the  funda¬ 
mental  solution  located  at  x  =  t  to  be  used  in  determining  the  poten¬ 
tial  of  the  nfth  Fourier  interference  component.  Let  gan*(0  be 
distribution  of  this  solution  along  the  body  axis  from  x  =  -1  to  x  =  co 
to  satisfy  the  boundary  conditions.  Then  the  potential  of  the  nTth 
Fourier  component  is 


cp2n(x,r,e)  =  J  cp2n*(x-£,r,e)  g2n*(t)  d! 

First,  the  set  of  fundamental  solutions  cp2n*(x,r,0)  will  be  presented, 
and  then  the  method  of  obtaining  gPp*(  5 )  will  be  discussed. 

Axial  strength  function.-  The  boundary  conditions  to  be  satisfied 
by  each  Fourier  Component of  the  interference  potential  are  given  by 
equations  (5)j  (6),  and  (9).  The  first  two  boundary  conditions  are 
satisfied  by  choosing  an  appropriate  set  of  fundamental  solutions.  Such 
a  set  of  fundamental  solutions,  derived  in  appendix  A,  is 


<P2n*(x-^r>0) 


cos  2nd  d211"1 
r2n  dx211"1 


[(x-| )2-r2] 


(11) 


The  fundamental  solution  is  located  at  the  point  x  =  £,  and  the  value 
of  £  varies  from  -1  to  +«o.  For  n  =  0  an  integration  is  to  be  per¬ 
formed.  The  last  boundary  condition  is  satisfied  by  distributing 
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fundamental  solutions,  as  given  in  equation  (ll),  along  the  body  axis  in 
variable  strength.  In  appendix  B  the  following  equation  is  derived  for 
determining  the  axial  strength  function 

g2n(x)  =  22n  rX  M2n(x-| )  f2n«  (5 )  d|  (12) 

( 4n i )J  o 

The  evaluation  of  gsn(x)  from  this  equation  is  done  numerically  or 
graphically.  The  values  of  f2nf(x)  are  obtained  on  the  basis  of  equa¬ 
tion  (10).  A  table  of  the  functions  M2n(x)  are  included  as  table  I  for 
use  in  equation  (12),  and  figure  5  illustrates  these  functions.  The 
properties  of  M2n(x)  are  discussed  in  reference  5» 

Potential  or  span  loading.-  Once  the  axial  strength  distribution  of 
the  fundamental  solution  is  known,  the  determination  of  the  interference 
potential  or  pressure  follows  readily.  Consider  the  general  point  P, 
shown  in  figure  6,  where  the  potential  is  to  be  determined.  The  funda¬ 
mental  solution  along  the  axis  from  -1  to  x-r  influences  the  potential 
at  P  which  is  given  as 


n  x-r  x2n- 1  _  1 

cp  (x,r,9)  =  f  |—  [(x- Sivrs  Ban.(l)  4S 

211  J_  1  dx211-1 

(13) 

By  making  a  unit  translation  of  the  variable  of  integration  so  that  the 
lower  limit  is  zero,  there  is  easily  obtained 


<P2n(x,r,e) 


cos  2nd  PX  r+1 

r2n  J0 


aan-i 

dx211"1 


[ (x-|+l)2-r2]2n”^  g2n(M  d|  (Ik) 

Equation  (l^)  serves  to  determine  the  potential  anywhere  in  the  flow 
field  once  the  axial  strength  function  g2n(x)  has  been  obtained.  Since 
the  integrated  lift  per  unit  span  up  to  a  given  streamwise  position  x 
is  proportional  to  cp  at  that  position,  equation  (ik)  also  gives  the 
span  loading. 


Pressure  coefficient.-  On  the  basis  of  linear  theory,  the  pressure 
coefficient  is 


2u  =  _2  dcp/dx 
V  ~  V 


(15) 


\ 


* 


Ik 
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With  the  aid  of  equation  (15)  it  is 

x-r+i 


-  ***2 


cos  2n 6 
V 


readily  shown  from  equation  (l^)  that 
^2n(T,r )  ga(x-r+l-T)  dT  (l6) 


where- 


K2n(T>r ) 


1  a2D[T(T+2r)]2n~l 

r2n  ^T2n 


The  procedure  in  obtaining  the  pressure  coefficient  is  first  to 
obtain  the  axial  strength  function  from  equation  (12)  by  numerical  or 
graphical  integration,  and  then  to  obtain  P2n  for  any  point  in  the 
flow  field  from  equation  (l6)  by  the  same  means. 

Body  pressure  coefficient.-  The  determination  of  the  pressure  coef¬ 
ficient  by  the  method  outlined  requires  two  numerical  or  graphical  inte¬ 
grations.  In  appendix  C  a  simpler  method  involving  only  one  integration 
is  derived  for  the  body  pressure  coefficient.  The  result  for  the  pres¬ 
sure  coefficient  is 


P2n  =  f  cos  2n0  f2n(x)  ~  |  cos  2n0  J  f2nU)W2n(x-£)  d5  (U) 

at  r  =  1.  The  functions  W2n(x)  are  universal  functions.  They  are 
tabulated  in  table  II,  and  are  plotted  in  figure  7* 

Combination  potential.-  The  combination  potential  is  obtained  by 
simply  adding  the  wing-alone  potential  to  the  interference  potential  cp. 
The  interference  potential  cp  is  obtained  as  already  described  by  add¬ 
ing  together  the  interference  potentials  of  a  number  of  Fourier  compo¬ 
nent  combinations.  The  lift  or  pressure  coefficients  for  the  combina¬ 
tion  are  similarly  obtained  simply  by  adding  to  those  for  the  wing  alone 
those  for  the  interference  combination. 

Referring  again  to  figure  2,  the  detailed  method  for  solving  either 
cases  (b)  or  (c)  has  been  summarized.  To  solve  the  general  case  of  both 
body  and  wing  at  angle  of  attack,  both  (b)  and  (c)  must  be  solved.  The 
wing  alone  for  each  case  is  different.  The  simplest  interference  problem 
is  that  for  the  body  at  zero  angle  of  attack.  Under  these  circumstances, 
is  zero,  and  only  case  (b)  must  be  solved.  This  case  will  be  calcu¬ 
lated  for  a  rectangular  wing -body  combination  as  an  illustrative  example. 
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PRESSURE  FIELD  ACTING  ON  A  COMBINATION  WITH  RECTANGULAR  WINGS 


One  of  the  simplest  wing-body  combinations  of  technological  impor¬ 
tance  is  that  formed  by  a  rectangular  wing  mounted  at  incidence  on  a  cir¬ 
cular  body  at  zero  angle  of  attack.  The  complete  pressure  field  acting 
on  such  a  combination  will  now  be  determined.  The  results  are  inter¬ 
preted  in  terms  of  lift  and  center  of  pressure  presented  as  a  function 
of  effective  aspect  ratio  and  effective  chord-radius  ratio.  As  pre¬ 
viously  mentioned,  the  wing  alone  can  be  specified  in  any  convenient 
manner  and,  for  the  purpose  of  the  example,  the  wing  alone  is  taken  as 
the  rectangular  wing  extending  straight  through  the  body  from  side  to 
side.  Although  the  analysis  as  carried  out  is  for  M  =  </2^  the  results 
are  presented  in  a  form  applicable  to  a  range  of  Mach  numbers.  The  steps 
in  performing  the  calculation  are:  (l)  to  determine  cp^,  the  wing-alone 
potential;  (2)  to  determine  the  velocity  amplitude  functions,  f^x); 

(3)  to  determine  the  axial  strength  functions,  g^x);  and  (4)  to  deter¬ 
mine  the  potential  or  pressure,  as  desired,  anywhere  in  the  field.  No 
tip  effects  are  considered  until  the  results  are  presented  as  a  function 
of  wing  aspect  ratio. 


Wing -Alone  Potential 


The  wing-alone  flow,  exclusive  of  tip  effects,  can  be  determined 
from  the  Ackeret  theory.  The  flow  at  a  spanwise  station  out  of  the 
region  of  influence  of  the  wing  tips  is  illustrated  in  figure  8.  The 
potential  for  the  flow  above  the  wing  is 


=  Vx  when  z>x 


(18) 


is  uniformly  -a^V 


=  Vx  +  a^V(x-z)  when  z<x  (19) 

— -  is  zero  and  the  downwash 
8y 

The  downwash  causes  a  flow  normal  to  the  sur- 


The  sidewash  produced  by  such  a  potential 
Sz 

face  r  =  1  in  amount  -o^V  sin  9.  This  means  that  for  a  body  conform¬ 
ing  to  the  wing-alone  flow,  the  deformation  is  zero  at  the  wing-body 
junctures  and  a  maximum  on  the  top  of  the  body.  The  interference  combi¬ 
nations  when  added  to  the  deformed  body  straighten  it  out  Fourier  com¬ 
ponent  by  Fourier  component. 


Fourier  Amplitudes  of  Body  Normal  Velocity 


The  Fourier  amplitudes  of  the  normal  velocity  induced  by  the  wing- 

89  u 

alone  potential  at  the  body  is  determined  by  expanding  — - 

Sr 


at  r  =  1 
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in  a  Fourier  cosine  series  of  even  multiples  of  0.  The  normal  velocity 
distribution  is  shown  in  figure  9.  For  x  >1  the  body  is  totally 
immersed  in  the  wing  downwash  field.  With  the  usual  equation  for  obtain¬ 
ing  the  Fourier  amplitudes  of  a  function,  there  is  obtained 


fo(x)  =  r 


sin_;Lx 


ctyV  sin  0  d0 


(20) 


fpn  (x )  -  ~ 


^  ^sin  ±x 


OyV  sin  0  cos  2n0  d0 


(21) 


The  integrations  give 


f0(x)  =  — “  (l  -  •/ L-x2 ^  when  x<l 
f0(x)  =  when  x>l 


(22) 

(23) 


Wx>  -  ^ 


2  cos  (2n-l)5  _  2  cos  (2n+l)5  4 

2n-l  2n+l  4n2-l  j 


^  /  \  -4Vaw 

^2n(x)  ^ — r  when  x^l 


when  x  <  1 


Jt(4n  -l) 


(24) 

(25) 


-1 


where  6  =  sin  x.  The  f2n(x)  functions  are  shown  in  figure  10.  The 
constant  values  of  f2n(x)  f°r  X~1  are  noteworthy. 


Axial  Strength  Functions 


The  calculation  of  the  axial  strength  distributions  has  been  carried 
out  with  the  aid  of  equation  (12).  The  integrations  were  performed 
numerically  by  means  of  integration  formulas  given  in  reference  12.  Two 
separate  cases  occur,  x<l  and  x>l.  For  x<l  a  typical  plot  of  the 
integrand  and  its  components  is  shown  in  figure  11(a).  It  will  be 
observed  that  the  Mq(x-| )  function  has  a  singularity  at  the  upper  limit 
of  integration.  Since  the  singularity  is  of  the  square  root  type,  it  is 
easily  handled.  For  x>l  the  determination  of  g0(x)  is  complicated 
by  the  fact  that  f0f(£)  has  a  square  root  singularity  at  1=1,  but 
simplified  by  the  fact  that  f0'(£)  is  zero  for  £  >1.  A  typical  plot 
of  the  integrand  and  its  components  for  x>l  is  shown  in  figure  11(b). 
The  singularity  of  MQ(x-l)  has  been  replaced  by  that  of  f0'(i)  at  the 
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upper  limit  of  integration.  For  x  =  1  the  square  root  singularities 
of  Mq(x-£  )  and  f0r(£)  both  occur  at  5  =  1>  reinforcing  each  other  in 
such  a  manner  that  gQ(x)  has  a  logarithmic  singularity  at  x  =  1.  In 
fact,  it  can  be  shown  that 

S2n(x)^^i 1  (-!)»  (i)  2S  log|l-x| 

(4ni)  \ny  k 

in  the  neighborhood  of  x  =  1.  The  function  ga(x)  minus  the  logarith¬ 
mic  term  passes  smoothly  through  x  =  1. 

A  plot  of  the  g3n(x)  functions  for  n  =  0  to  n  =  3  is  shown  in 
figure  12.  The  main  effect  noted  is  that  the  number  of  times  the  curve 
crosses  the  zero  axis  increases  as  the  order  of  the  g^Cx)  functions 
increases.  The  singularities  at  x  =  1  are  also  illustrated. 


Pressure  Distributions 


Obtaining  the  strength  functions  of  g2Q(x)  is  tantamount  to  getting 
the  potential  or  pressure  coefficient  anywhere  in  the  flow  field  as  given 
by  equation  (l4)  or  (l6).  If  only  the  span  load  distribution  is  required, 
then  one  calculation  using  equation  ( 14 )  suffices  for  each  spanwise 
station;  but  if  the  detailed  pressure  distribution  is  required,  then 
many  calculations  using  equation  (l6)  must  be  made  for  each  station.  The 
detailed  pressure  distributions  have  been  calculated  using  equations  (l6) 
and  (17).  In  equation  (16),  the  factor 


(T(r+2r)]2n"i 


appears  inside  the  integral  as  an  influence  coefficient  for  converting 
axial  strength  function  into  pressure  coefficient.  The  factor  plays  an 
important  role  in  determining  the  accuracy  of  the  numerical  work.  In 
figure  13  plots  of  KQ(r,r)  and  K2(i-,r)  are  presented  showing  the  singu¬ 
larities  of  K2ri(T,r)  at  the  origin  of  t.  These  singularities,  together 
with  those  of  ga(x-r+l-T)  at  t  =  x-r,  bring  about  two  different  cases 
in  the  determination  of  the  pressure  coefficient.  If  x<r,  only 
the  K2n(T,r)  singularity  occurs  in  the  integration  as  illustrated  by 
figure  14(a).  However,  for  x>r  there  are  two  singularities  as  shown 
in  figure  14(b).  For  x  =  r  the  two  singularities  come  into  confluence, 
but  they  still  produce  a  finite  pressure  coefficient. 

Interference  pressure  distribution.-  The  interference  pressure  dis¬ 
tributions  have  been  calculated  for  the  first  four  Fourier  components 
and  are  presented  in  figure  15.  In  this  figure  the  abscissa  is  propor¬ 
tional  to  distance  behind  the  Mach  line  originating  at  the  leading  edge 
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of  the  juncture,  as  illustrated  in  the  figure.  Although  the  calcula¬ 
tions  have  been  carried  out  for  M  =  «/2,  that  is,  0=1,  and  for  unit 
radius,  they  are  generalized  to  all  Mach  numbers  and  body  radii  by  replac¬ 
ing  x-r+1  by  — - -  +  1  and  P2n  by  0P2n  as  has  been  done  in  the 

pa  a  ~  A 

figure.  From  the  figure  it  is  apparent  that  the  cusps  in  the  pressure 
distributions  are  propagated  downstream  along  lines  of  constant 
- -  —  +  1  or  x-0r;  that  is,  along  the  downstream  characteristics.  As 

the  pressure  distributions  move  outward  from  the  body  along  the  down¬ 
stream  characteristics,  they  are  distorted  and  decreased  in  magnitude. 


Increasing  the  order  of  the  Fourier  harmonics  causes  two  important 
effects:  first,  the  number  of  points  of  zero  pressure  is  increased  and, 
second,  the  pressure  coefficient  damps  more  rapidly.  As  a  result  of  the 
first  effect  the  contributions  of  the  higher  harmonics  to  the  combination 
span  loading  are  proportionally  less  than  their  contributions  to  the 
pressure  coefficient;  while  as  a  result  of  the  second  effect  the  more 
remote  a  point  is  from  the  leading  edge  of  the  wing -body  juncture,  the 
fewer  the  number  of  Fourier  components  that  must  be  included  to  obtain 
its  pressure  coefficient  accurately-  All  interference  pressure  distribu¬ 
tions  exhibit  discontinuities  in  slope  at  -£ - £  +  1  =  1.  This  behavior 

Ba  a 

is  a  consequence  of  the  fact  that  the  body  becomes  totally  immersed  in 
the  wing-alone  flow  field  for  this  condition.  When  the  pressure  distri¬ 
butions  of  the  various  Fourier  components  are  added  together  to  obtain 
the  interference  pressure  distributions,  the  discontinuities  in  slope 
tend  to  cancel  so  that  the  pressure  distribution  for  the  combination  will 
be  smooth. 


A  detailed  examination  of  the  interference  pressure  distribution 
for  the  first  Fourier  component  illustrates  several  points  of  interest. 
The  importance  of  the  component  arises  from  the  fact  that  it  accounts 
for  most  of  the  effect  of  interference  on  the  span  loading.  The  reason 
for  this  is  that  the  pressure  coefficients  for  n  =  0  are  of  invariable 
sign.  The  effect  of  the  first  Fourier  component  is  to  reduce  the  veloc¬ 
ity  induced  normal  to  the  body  by  the  wing-alone  flow  field  to  zero 
average  around  the  body  for  6  =  0  to  6  =  n  at  any  streamwise  location. 


For  purposes  of  comparison  with  the  exact  results  for  n  =  0,  some 
approximate  results  have  been  included  in  figure  15(a)*  For  values 

body,  the  Acker et  value  of  PQ  (twice  the  local  stream 
(3)  is  a  close  approximation  to  the  true  pressure  coef- 


of  — <1  on  the 

3a 

angle  divided  by 

ficient.  This  is  the  result  of  the  facts  that  the  part  of  the  body 
affecting  the  interference  is  effectively  plane  for  points  near  the  lead¬ 
ing  edge  of  the  wing -body  juncture  and  that  there  is  no  variation  of  any 
quantities  with  6  so  that  an  approximate  two-dimensional  situation  pre¬ 


vails  .  As 


3a 


increases  beyond  unity  on  the  body,  there  is  a  rapid 


decrease  in  the  pressure  coefficient  below  the  Ackeret  value  due  to  the 
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effect  of  all  disturbances  in  front  of  the  point  in  question  as  repre 
sented  by  the  integral  of  equation  (17). 


In  reference  5,  the  following  approximate  results  were  obtained  for 


small  and  large  values  of 


Gs-i  +  1) 


for  the  pressure  coefficient: 


0P 


2n 


cos  2n6 
*-/r~ 


— —  —  +  1 
Ba  a 


(26) 


BP  — J2L.  JL 
0  Jtx/Ba*  Ba 


£  +  l->  » 
a 


(27) 


Figure  15(a)  shows  that  for  values  of  ^  ~  f  +  1  <0.6  equation  (26)  is 
a  good  approximation  for  n  =  0  although  it  is  of  little  value  for 
higher- order  harmonics.  There  is  a  general  tendency  of  PQ  to  approach  a 
uniform  value  independent  of  r  as  ~  +  1  becomes  large,  as  shown 

by  equation  (27).  The  damping  in  the  characteristic  direction,  although 
initially  inversely  proportional  to  the  square  root  of  r,  is  ultimately 
independent  of  r. 

The  pressure  distributions  of  figure  15  were  determined  by  equa¬ 
tion  (17)  for  r  =  1,  and  by  equations  (12)  and  (16)  for  r  =  1.5,  2.0, 
and  3.0.  However,  the  values  for  r  =  1  were  checked  by  equations  (12) 
and  (16)  before  these  equations  were  applied  to  the  higher  values  of  r. 
For  r  =  1  equation  (17)  gives  more  accurate  results  than  equations  (12) 
and  (16)  since  it  requires  only  one  numerical  integration.  The  accuracy 
of  both  methods  of  calculation  decreases  as  n  and  —  —  +  1  increase 
because  the  small  values  of  the  pressure  coefficients  under  these  condi¬ 
tions  are  the  result  of  large  counterbalancing  influences . 

Pressure  distribution  in  juncture  of  wing-body  combination. -  By 
adding  the  interference  pressure  coefficients  of  the  various  Fourier 
components  to  that  for  the  wing  alone,  the  pressure  distribution  for  the 
combination  is  obtained.  The  addition  has  been  carried  out  for  the  wing— 
body  juncture  using  four  Fourier  components  and  six  Fourier  components, 
and  the  results  are  presented  in  figure  16.  The  pressure  coefficient 
with  interference  is  less  in  magnitude  than  two,  the  value  without  inter¬ 
ference,  showing  that  significant  losses  of  lift  occur  in  the  wing-body 
juncture.  A  comparison  of  the  results  for  four  components  and  six  com¬ 
ponents  shows  that  four  components  give  good  over-all  accuracy  for  all 
values  of  —  greater  than  1.  For  small  values  of  -jj~  in  the  wing- 

body  juncture,  the  curvature  of  the  body  insofar  as  the  flow  is  concerned 
is  not  large  so  that  the  body  is  effectively  a  vertical  boundary  on  which 
a  given  distribution  of  normal  velocity  is  producing  an  interference 
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field.  Supersonic  wing  theory  applied  to  this  condition  gives  for  the 
net  interference  pressure  coefficient  (reference  5) 

££  when  •£.  0  (28) 

ofy  3*pa  pa 

It  is  clear  that  the  calculated  results  can  he  joined  smoothly  to  this 
result-  Using  the  result  of  equation  (28)  enables  satisfactory  results 
to  be  obtained  with  four  Fourier  components. 

The  critical  region  in  the  convergence  of  the  solution  is  that  near 

the  leading  edge  of  the  wing-body  juncture.  The  higher  harmonics  have 

their  most  important  effect  near  here ,  and  rapidly  damp  downstream  along 

the  body.  Hence  more  and  more  Fourier  components  would  be  required  to 

get  accuracy  for  smaller  and  smaller  values  of  However,  with  the 

pa 

result  of  equation  (28),  this  extra  work  is  unnecessary. 

One  point  of  interest  in  figure  l6  is  the  fact  that  when  equals 

approximately  3>  the  pressure  coefficient  increases  in  magnitude.  This 

is  due  to  the  fact  that  for  — >*  the  influence  of  the  opposite  half- 

Pa 

wing  is  felt  in  the  wing-body  juncture. 

Pressure  distribution  on  top  meridian  of  wing-body  combination. - 
The  pressure  distribution  on  the  top  meridian  of  the  wing-body  combina- 
tion  is  obtained  in  the  same  fashion  as  that  at  the  wing -body  juncture, 
the  difference  being  that  the  pressures  due  to  the  even  number  Fourier 
components  have  the  same  sign  at  the  meridian  as  at  the  juncture,  whereas 
the  odd  numbered  components  have  reversed  signs.  The  pressure  distribu¬ 
tions  based  on  four  and  six  Fourier  components  are  shown  in  figure  17 • 

Several  interesting  effects  are  exhibited  by  the  results.  The  step 

in  the  wing-alone  pressure  at  =  1  is  effectively  canceled  by  the 

interference  pressures  of  the  Fourier  components  from  -j~  =  1  to  =  */£> 

and  for  —  >it/2  the  pressure  increases  rapidly  and  tends  toward  the  two- 
Pa 

dimensional  value.  The  effect  of  the  interference  pressure  in  canceling 

the  effect  of  the  wing  alone  on  the  top  of  the  body  from  =  1 

to  —  =  rt/2  is  to  be  expected  since  the  wing  of  the  combination  can 
pa 

have  no  effect  on  the  top  of  the  body  unless  as  has  been 

pa 

already  pointed  out.  If  an  infinite  number  of  Fourier  components  had 
been  taken,  the  pressure  coefficients  would  be  identically  zero 

from  —  =  0  to  —  =  it/2.  The  general  behavior  in  this  regard  is  evi- 
pa  Pa  ' 

dence  of  the  plausibility  of  the  calculated  results. 
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The  tendency  of  the  pressures  to  approach  an  asymptotic  value  is  also 
illustrated  by  figure  17 •  This  asymptotic  value  represented  by  the  sum 
of  the  ving-alone  pressure  plus  the  asymptotic  results  for  the  first 
Fourier  component  is  given  by  the  following  equation: 


P*-2  + 


k 

flx/pa 


when 


>2.4 


For  —  >2»k,  the  results  of  this  equation  are  in  good  agreement 
with  the  results  of  figure  17* 

Some  evidence  is  furnished  from  the  pressure  calculations  for  the 
juncture  and  top  of  the  body  concerning  the  number  of  Fourier  components 
necessary  for  accuracy.  Comparisons  made  in  figures  1 6  and  17  show  that 
about  four  components  are  sufficient,  and  that  the  addition  of  two  more 
is  not  worth  the  extra  work. 

Pressure  distributions  for  wing  of  wing-body  combination.-  The  dis¬ 
tribution  of  the  pressure  acting  on  the  wing  of  the  combination  can  be 
determined  in  a  manner  similar  to  that  for  the  wing-body  juncture  by 
adding  to  the  wing-alone  pressure  those  due  to  the  Fourier  components. 

The  resultant  pressure  distribution  for  the  wing  based  on  four  Fourier 
components  is  shown  in  figure  l8.  A  drawing  of  the  pressure  field  for 
the  complete  configuration  is  included  as  figure  19*  For  small  values 
of  — -  the  higher -order  oscillations  in  the  pressure  coefficient  as 

Del 

shown  in  figure  1 6  have  been  ignored  and  the  curves  have  been  faired 
through  them. 

Since  the  region  of  influence  of  the  body  on  the  wing  is  confined  to 
the  wing  region  downstream  of  the  Mach  lines  emanating  from  the  leading 
edge  of  the  juncture,  in  front  of  this  line  the  pressures  are  uniform 
at  the  two-dimensional  value  and  behind  the  line  there  is  a  decrease  in 
the  magnitude  of  the  pressure  coefficient.  If  the  body  were  a  perfect 
reflector,  that  is,  a  vertical  wall  of  infinite  extent,  then  there  would 
be  no  pressure  loss.  However,  the  pressure  pulses  originating  on  the 
wing  are  only  in  part  reflected  by  the  circular  body.  The  efficiency  of 
the  body  as  a  reflector  is  discussed  subsequently  in  connection  with  span 
loading.  The  tendency  of  the  pressure  to  increase  in  magnitude  near  the 
inboard  trailing  edge  is  due  to  the  effect  of  the  opposite  wing  panel 
which  at  the  wing  juncture  is  felt  downstream  of  the  point  =  jt. 


Span  Loading 


The  span  load  distributions  for  a  range  of  rectangular  wing -body 
combinations  with  the  body  at  zero  angle  of  attack  can  be  determined 
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from  the  pressure  distributions  of  figure  18.  For  purposes  of  specify¬ 
ing  the  span  loading  the  following  equation  is  used; 

■\ 


L 


!a*)ay 


(30) 


The  span  loading  is  taken  as  the  quantity  inside  the  brackets. 


The  pressure  results  of  figure  18  are  for  values  of  the  effective 

chord-radius  ratio  of  4  or  less  and  for  values  of  the  effective  aspect 

ratio  of  2  or  greater.  Span  loadings  for  any  combination  of  —  (or  c*) 

^a 

and  pA  in  these  ranges  can  be  obtained  by  integrating  the  pressure  dis¬ 
tributions.  The  span  loading  evaluations  have  been  made  for  c*  =  4 
and  pA>2.  First  the  span  loadings  due  to  the  various  Fourier  components 
are  discussed,  and  then  the  span  loadings  for  the  actual  wing-body  combi¬ 
nations  are  presented. 


Fourier  component  combinations.-  In  figure  20,  the  contributions  to 
the  span  loading  for  the  first  three  Fourier  components  are  shown. 

For  n=  0  the  pressure  field  does  not  depend  on  0,  being  axially  sym¬ 
metric,  and  a  constant  loading  exists  on  the  body.  However,  on  the  wing 
as  the  spanwise  distance  increases  there  is  a  decrease  in  the  span  load¬ 
ing  due  primarily  to  decrease  in  the  length  of  chord  over  which  the 
interference  pressures  act.  The  span  loading  due  to  the  first  Fourier 
component  causes  a  loss  of  lift  everywhere  along  the  span. 

A  comparison  of  the  results  of  figure  20  for  n  =  0  and  n  =  1  shows 
that  the  first  Fourier  component  accounts  almost  entirely  for  the  effect 
of  interference  on  the  span  loading  of  the  combination.  For  the  body 
this  fact  is  even  more  true  than  for  the  wing.  This  fact  is  of  consider¬ 
able  importance  since  it  gives  a  simple  means  of  extending  the  lift  and 

moment  results  to  larger  values  of  than  those  for  which  the  pressure 

pa 

distributions  have  been  calculated.  Also,  it  suggests  a  simple  means  of 
minimizing  the  adverse  effects  of  interference  on  lift  as  will  subse¬ 
quently  be  pointed  out. 


With  the  techniques  of  Laplace  transform  theory,  it  is  possible  to 
obtain  asymptotic  formulas  for  the  span  loadings  of  the  various  Fourier 
components.  For  the  first  Fourier  component  the  following  asymptotic 
result  has  been  obtained  by  the  standard  methods  of  Laplace  transform 
theory.  (See  reference  5*) 


-2 


dx-»  |  log 


X 

Pa 


* 


when 


— >  00 


(31) 


V 


NACA  TN  2677 


23 


The  asymptotic  result  for  the  span  loading  given  by  this  equation,  when 
compared  with  the  results  of  the  exact  calculations  in  figure  20,  is 
seen  to  be  slightly  high.  However,  for  values  of  £-  greater  than  4 

pel 

the  difference  between  the  results  decreases,  and  equation  (31)  thus  pro¬ 
vides  a  satisfactory  means  of  extrapolating  the  results  of  the  present 

calculations  for  span  loading  to  larger  values  of  — . 

w  0a 


The  asymptotic  result  has  also  been  determined  for  the  higher-order 
Fourier  components  as  a  matter  of  interest.  The  span  loading  is 


Pan  ^  ^  8  cos  2n9 

“W  jtn^n2-!)!-211 


16  ^r2n  +  —4^)  (4n-l)lcos  2n0 
*(2nl  )2(4n2-l)24n“1  (^T 


as  -2«- 
0a 


(32) 


A  comparison  between  the  results  of  equation  (32)  and  the  exact  solution 
for  n  =  1  in  figure  20  shows  agreement  and  corroborates  the  fact  that 
the  span  loading  of  all  but  the  first  Fourier  component  is  negligible 
l°r  It  is  also  to  be  noted  that  the  contribution  to  the  loading 

of  the  first  component  given  by  equation  (31)  increases  without  limit 
as  x  — >  °°>  whereas  the  span  loading  of  the  higher-order  components  is 
finite . 


Family  of  complete  combinations.-  To  obtain  the  span  loading  for  the 

family  of  combinations  for  which  —  =  4,  it  is  necessary  to  consider  the 

pa 

loading  of  both  "the  wing  alone  and  the  Fourier  components  •  The  necessary 
calculations  have  been  carried  out  and  the  span  loadings  for  the  family 
of  combinations  based  on  one  and  four  Fourier  components  are  both  shown 
in  figure  21.  The  loading  due  to  the  wing  alone  is  also  shown.  No 
effect  of  wing  tips  has  been  included.  It  is  to  be  noted  in  figure  21 
that,  whereas  the  loading  on  the  wing  due  to  its  own  pressure  field  is 
constant,  there  is  some  loss  on  the  body  because  of  the  fact  that  the 
pressure  field  of  the  wing  alone  acts  on  the  body  only  if  x>pa  sin  0. 
However,  if  an  afterbody  is  included,  some  of  the  lift  lost  can  be  recov¬ 
ered.  As  has  already  been  pointed  out,  the  pressures  due  to  the  first 
Fourier  components  are  positive  on  the  upper  half  of  the  wing-body  combi¬ 
nation  and  produce  a  loss  of  lift,  as  figure  21  shows.  When  the  effects 
of  four  Fourier  components  are  taken  into  account,  the  net  lift  is 
slightly  higher  than  that  for  one  Fourier  component,  but  the  difference 
is  not  significant.  For  most  engineering  purposes,  one  Fourier  component 
is  sufficient  for  determining  the  span  loading  when  — >  h. 

pa 

Some  insight  into  the  mechanism  of  wing-body  Interference  can  be 
gained  by  comparing  the  span  loading  for  the  combination  with  those  for 
two  reference  loadings:  (l)  the  complete  reflection  case  for  which  the 
blanketed  area  of  the  wing  acts  effectively  at  ay,  and  (2)  the 
no-reflection  case  for  which  the  blanketed  area  of  the  wing  is  supposed 
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to  act  effectively  at  zero  angle  of  attack.  The  span  loading  correspond¬ 
ing  to  the  first  case  of  complete  reflection  of  the  wing  pressure  pulses 
by  the  body  is,  in  fact,  the  span  loading  marked  "wing  alone"  in  fig¬ 
ure  21.  A  comparison  of  this  curve  with  that  based  on  one  or  four 
Fourier  components  shows  that  the  loading  given  on  the  assumption  that 
the  wing  blanketed  area  is  fully  effective  in  lift  is  too  optimistic. 
Under  the  conditions  of  the  second  reference  loading,  the  sole  purpose  of 
the  blanketed  area  is  to  support  lift  generated  by  the  wings.  A  compari¬ 
son  of  the  span  loading  for  this  case  with  the  true  loading  shows  that 
the  average  load  on  the  body  is  well  predicted,  but  that  the  loading  on 
the  wing  is,  of  course,  underestimated.  A  comparison  of  the  true  load¬ 
ing  with  those  for  the  two  reference  cases  reveals  the  interesting  fact 
that  the  body  is  somewhat  less  than  50  percent  effective  in  reflection 
for  this  particular  family  of  configurations. 


Lift 


the  pressure  distributions  already  presented 


For  values  of  —  <4 

pa 

are  sufficient  for  obtaining  span  loading  or  lift  on  either  the  wing  or 

body  for  all  combinations  having  sufficiently  large  aspect  ratios  to 

avoid  effects  of  the  tips  on  the  wing -body  interference.  This  is  the 

case  for  pA>2.  For  values  of  larger  than  four,  the  lift  results 

pa 

are  presented  in  terms  of  a  nondimensional  parameter  kw,  defined  as  the 
ratio  of  the  lift  on  the  exposed  half -wings  in  combination  (exclusive  of 
that  on  the  body)  to  that  on  the  exposed  half -wings  joined  together. 


kw  = 


LWC 

Lw 


(33) 


For  ~>k  the  value  of  kw  can  be  obtained  by  using  the  asymptotic 
form  of  the  span  loading  given  by  equation  (31). 
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c*  — >  00  (34) 

The  values  of  kv  have  been  determined  from  the  pressure  distribu¬ 
tions  of  figure  18  for  values  of  JL<4  and  from  equation  (34)  for 

pa 

values  of  —  >4.  The  effect  of  the  wing  tips  has  been  taken  into 
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account  by  utilizing  reference  13.  The  results  are  shown  in  figure  22 
wherein  kw  is  given  as  a  function  of  ~  for  various  effective  aspect 

ratios  of  2  and  greater.  It  should  be  borne  in  mind  that  the  results  of 
the  figure  are  for  a  combination  of  body  and  rectangular  wing  or  an  all¬ 
movable,  rectangular  control  surface  with  no  gap.  It  is  noted  in  the 
figure  that  the  exact  results  for  can  ^e  Tailed  into  the  asymptotic 

results  for  thereby  providing  a  design  chart  for  engineering  pur- 

pa  c 

poses  for  the  entire  range  of  The  curves  of  figure  22  illustrate  the 

decrease  of  kw*  as  increases  at  constant  effective  aspect  ratio, 

and  the  slow  increase  of  kv  toward  unity  as  the  wing  chord  becomes 
very  large.  The  loss  of  lift  is  most  serious  for  pA  =  2,  being  about 
15  percent  in  the  worst  case. 

A  practical  point  in  connection  with  the  loss  of  lift  on  the  wing 
due  to  interference  is  that  this  loss  occurs  no  matter  what  the  body 
angle  of  attack,  even  though  the  calculations  are  made  for  cx^  =  0.  It 
occurs  either  in  the  case  of  a  wing  mounted  on  a  body  or  in  the  case  of 
a  deflected  all -movable  control  surface.  For  wings  with  swept  leading 
edges  for  which  all  of  the  wing  area  lies  in  the  region  affected  by  the 
interference,  even  larger  losses  than  occur  with  rectangular  wings  are  to 
be  anticipated.  However,  the  loss  of  lift  at  the  design  condition  can, 
at  least  in  principle,  be  largely  prevented  by  designing  the  fuselage 
so  that  it  conforms  to  the  first  Fourier  component  in  the  wing-alone 
flow.  This  would  involve  contracting  the  fuselage  above  the  horizontal 
plane  of  symmetry  in  a  rotationally  symmetric  fashion  and  expanding  a 
like  amount  beneath  the  horizontal  plane  of  symmetry.  Whether  or  not  . 
such  a  change  would  improve  the  lift-drag  ratio  can  best  be  determined 
by  experiment. 


Center  of  Pressure 


The  center-of -pressure  locations  have  been  calculated  for  the  same 
condition  as  the  lift  results  of  figure  22.  The  center-of -pres sure  loca¬ 
tion  in  chord  lengths  behind  the  leading  edge  are  presented  in  figure  23* 
For  large  values  of  c*,  an  asymptotic  result  has  been  calculated  for 
xc.p./c  usinS  the  methods  of  Laplace  transform  theory  and  considering 
only  one  Fourier  component. 


1 

3PA 


J-Qg  2c*  +  JL 

2c*2  jtc* 


1 

(3  Ac*3 


as  c*  — >  oo 


(35) 


2  6 
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The  values  of 


xc.p. 


have  been  determined  from  the  pressure  distri¬ 


butions  of  figure  18  for  values  of 


and  by  equation  (35)  for 

values  of  —->4.  The  loss  of  lift  near  the  tips  has  been  taken  into 


pa 

consideration. 


The  exact  results  for 


-£-<4 

pa 


have  been  faired  into  the 


asymptotic  results  for  large  values  of  by  dashed  curves  to  provide 


Pa 


It  is 


corre- 


an  engineering  design  chart  covering  the  entire  range  of 

Pa 

again  mentioned  that  this  chart  is  applicable  both  to  the  wing  of  an 
airplane  or  missile  or  to  an  all-movable,  rectangular  control  surface 

with  no  gap.  The  curves  of  figure  23  start  at  values  of  ~ *■£ 

c 

sponding  to  those  for  the  wing  alone  at  -iL  =  0.  As  JL  increases  for 

Pa  Pa 

constant  pA  there  is  a  forward  movement  of  the  center  of  pressure 
because  of  the  loss  of  lift  due  to  interference  which  is  mostly  effective 
on  the  rear  of  the  wing.  For  the  lowest  effective  aspect  ratio  of  2 
there  is  about  a  4-percent  forward  movement  of  the  center  of  pressure 
due  to  interference  in  the  extreme  case.  For  large  effective  aspect 

ratios  the  forward  movement  is  not  nearly  so  large.  As  the  value  of 

□cl 

increases  for  constant  pA  there  is  an  asymptotic  approach  of  the  center 
of  pressure  back  to  the  wing-alone  value. 


Drag 


The  results  of  the  lift  calculations  can  be  used  directly  to  deter¬ 
mine  the  effect  of  interference  on  the  wave  drag  at  zero  angle  of  attack 
of  a  symmetrical,  double-wedge,  rectangular  wing  mounted  on  a  circular 
cylinder.  In  fact,  it  is  not  difficult  to  derive  the  result  that  this 
drag  coefficient  based  on  the  area  of  the  exposed  wing  is  given  simply  as 


PCp 

(t/c)2 


=  8 


1 


MpA,c*) 


(36) 


Equation  (38)  is  for  the  maximum- thickness  position  at  the  midchord, 
although  other  maximum-thickness  positions  can  easily  be  handled. 


The  values  of  the  drag  parameter  given  by  equation  (36)  have  been 
calculated  for  effective  aspect  ratios  of  2,  3,  and  4  using  the  values 
of  k^  from  figure  22.  The  results  of  the  calculations  are  shown  in 

P°D  c 

- -  is  plotted  against  —  for  constant  values 

(t/c  f  pa 

_  0  PCD 
Pa  (t/c)‘ 


figure  2k  wherein 


of  PA.  At 


is  equal  to  the  wing-alone  value  of  4. 
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As  ■£-  increases,  there  is  an  increase  in  drag  which  in  the  extreme  case 
pa 

amounts  to  about  4  percent  for  J3A  =  2,  and  thereafter  there  is  a  steady 

decrease  in  drag  until  at  -p-  =  12,  the  upper  range  of  the  calculations , 

pa 

there  is  about  a  4-percent  decrease  in  drag  from  the  wing-alone  value. 

For  other  effective  aspect  ratios  the  changes  are  inversely  proportional 
to  the  effective  aspect  ratio.  It  can  thus  be  said  for  rectangular 
wings  of  the  present  type  with  |3A>2  and  —  <12  that  no  important 

pel 

effects  of  interference  on  the  minimum  drag  occur.  For  missile-like  con¬ 
figurations  for  which  the  wing  contributes  only  a  small  part  of  the  drag, 
the  effects  of  interference  on  minimum  drag  will  be  negligible.  For 
wing-body  combinations  having  a  large  part  of  the  wing  lying  behind  the 
Mach  waves  from  the  leading  edge  of  the  wing-body  junctures,  there  can 
be  appreciable  effects  of  interference  on  minimum  wave  drag. 


CONCLUDING  REMARKS 


A  theoretical  method  has  been  developed  that  makes  it  possible  to 
obtain  the  pressure  field  acting  on  a  wing-body  combination  composed  of 
a  wing  with  supersonic  leading  and  trailing  edges  and  a  circular  body. 

The  method  has  been  applied  to  the  calculation  of  the  pressure  fields 
acting  on  combinations  of  bodies  and  rectangular  wings  with  the  bodies  at 
zero  angle  of  attack.  The  exact  calculations  are  for  combinations  for 
which  the  effective  aspect  ratio  of  the  exposed  wing  panels  joined  together 
is  greater  than  2,  and  for  which  the  effective  chord-radius  ratio  is  4  or 
less.  For  the  family  of  combinations  for  which  the  effective  chord-radius 
ratio  is  4,  it  was  found  that  the  blanketed  area  of  the  wing  acts  effec¬ 
tively  at  zero  angle  of  attack  and  served  primarily  to  support  lift  carried 
over  from  the  wing  onto  the  body.  It  was  determined,  that  for  this  family 
of  combinations  the  body  was  somewhat  less  than  50  percent  effective,  com¬ 
pared  to  a  perfect  reflection  plane,  in  reflecting  pressure  waves  back 
onto  the  wing. 

By  comparing  the  contributions  to  the  span  loading  of  the  various 
Fourier  components,  it  was  determined  that  for  rectangular  wing -body 
combinations  with  effective  chord-radius  ratios  of  4  or  greater,  the 
loss  of  lift  due  to  interference  is  associated  principally  with  the  first 
component.  On  the  basis  of  this  fact  it  is  possible  to  obtain  simple 
asymptotic  formulas  for  wing  lift  and  wing  center -of -pres sure  positions 
for  large  values  of  the  effective  chord-radius  ratio.  The  results  of 
the  exact  calculations  and  of  those  based  on  the  asymptotic  formulas 
enabled  the  construction  of  design  charts  showing  the  wing  lift  and 
center -of' -pressure  position  as  a  function  of  effective  aspect  ratio  and 
effective  chord-radius  ratio.  These  charts  are  applicable  to  wings  or 
control  surfaces  deflected  in  the  presence  of  a  round  body.  The  charts 
show  that  as  a  result  of  interference  the  lift  of  the  wing  panels  in  com¬ 
bination  can  be  reduced  as  much  as  15  percent  from  that  of  the  half -wings 
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joined  together ,  and  that  the  center  of  pressure  of  the  wing  can  move 
forward  as  much  as  k  percent  of  the  chord.  No  large  effects  of  inter¬ 
ference  on  the  wave  drag  of  symmetrical  double-wedge  rectangular  wings 
with  the  maximum  thickness  at  the  midchord  were  found. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  January  8,  1952. 
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* 

appendix  a 

DETERMINATION  OF  FUNDAMENTAL  SOLUTIONS  FOR  FOURIER 
COMPONENT  COMBINATIONS 


Von  Karman  and  Moore  determine  the  potential  for  a  body  of  revolu¬ 
tion  by  distributing  along  the  axis  of  the  body  fundamental  solutions  of 
an  axially  symmetric  type.  This  naturally  suggests  the  use  of  fundamen¬ 
tal  solutions  varying  as  cos  2n0  for  the  Fourier  component  combinations 
for  which  the  potential  varies  as  cos  2n0.  Such  fundamental  solutions 
satisfy  the  boundary  condition  given  by  equation  (5)*  The  boundary  con¬ 
dition'  given  by  equation  (6)  can  be  satisfied  by  making  the  fundamental 
solution  imaginary  in  the  forecone,  that  is,  by  making  it  proportional 
to  [(x-£ )2-r2]m/2  where  m  is  an  odd  integer.  These  conditions  make 
the  desired  fundamental  solutions  of  the  form 

cos  2n0  [  (x-|  )2-r2]m//2  ^(x-|,r) 

where  t(x-5, r)  Is  a  function  to  be  determined  in  such  a  manner  that  the 
fundamental  solution  obeys  the  wave  equation.  The  parameter  £  is  the 
position  on  the  x  axis  where  the  fundamental  solution  is  located.  In 
the  case  of  a  wing-body  combination,  £  will  vary  from  -1  to  +00.  How¬ 
ever,  a  value  of  5=0  is  used  in  the  derivation  of  the  fundamental 
solution  so  that  It  will  be  at  the  origin.  If,  as  a  simple  possibility, 
i)r(x,r)  is  taken  as  a  power  of  r  (or  x)  and  the  fundamental  solution  is 
substituted  into  the  wave  equation  in  cylindrical  coordinates  (see  fig.  1 
for  coordinate  system) 


cos  2n0 


cos  2n0 


(Al) 

solutions  are  obtained 

(x2-r2)2n~i 

r2n 

(A2) 

r2n 

- ~r 

(A3) 

(x2-r2)2n+2 


The  first  solution  is  singular  along  the  body  axis  (r  =  0),  and  the 
second  solution  is  singular  along  the  Mach  cone  (x  =  r).  Derivatives  of 
the  above  solutions  with  respect  to  x  are  again  possible  fundamental 
solutions. 

A  set  of  fundamental  solutions  having  its  singularities  along  the 
body  axis  would  be  better  than  a  set  having  its  singularities  on  the 
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Mach,  cone  since  it  was  desired  to  get  the  external  flow  field  about  a 
combination.  However,  when  the  foregoing  fundamental  solutions  with  the 
singularities  on  the  axis  were  used,  it  was  discovered  that  the  functions 
defining  the  strengths  of  the  distributions  along  the  body  axis  had 
unmanageable  singularities.  In  fact,  only  by  using  the  (2n-l),th  deriv¬ 
ative  with  respect  to  x  of  these  solutions  was  another  set  of  funda¬ 
mental  solutions  obtained  which  gave  manageable  singularities  in  the 
strengths  of  the  distributions  along  the  axis.  The  set  of  fundamental 
solutions,  cp2n*(x,r,e), 

.  £2^22  gi.  (x=-r (M) 

are  those  adopted  for  this  report.  For  n  =  0,  the  partial  differentia¬ 
tion  convention  of  equation  (Ak)  is  taken  to  indicate  an  integration 
with  regard  to  x. 
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APPENDIX  B 

DETERMIMTI ON  OF  AXIAL  STRENGTH  FUNCTIONS  FOR  FOURIER 
COMPONENT  COMBIMTIONS 


The  distribution  along  the  axis  of  the  fundamental  solution  Is  deter¬ 
mined  by  the  boundary  condition  of  equation  (9),  that  is,  the  distribu¬ 
tion  is  to  produce  a  given  velocity  amplitude  function  f2n(x).  If  the 
distribution  of  cp2n*(x,r,0)  (equation  (A4) )  along  the  body  axis 
from  -1  to  <*>  to  produce  a  step  function  of  f2n(x)  were  known,  then 
the  axial  strength  distribution  to  produce  the  given  function  f2n(x) 
could  be  closely  approximated  by  step  functions  as  shown  in  figure  25. 

Let  h2n(x-l)  be  the  distribution  along  the  axis  from  -1  to  a  to  pro¬ 
duce  a  unit  step  in  f2n(x)  at  the  origin.  Then  to  obtain  the  axial 
strength  function  to  produce  the  stepwise  approximation  to  f2n(x)  shown 
in  the  figure,  the  following  sum  must  be  formed! 


&3+1 


If  g2n*(x-l)  is  the  strength  distribution  along  the  axis  for  the  actual 
f2n(x)  distribution  found  by  letting  the  Ax's  approach  zero,  then 
since  ■^2Xl(0)  "  0 


g2n*(x-l)  =  J  Wx-l-t)f2n'(0  d.§ 


(Bl) 


vith  the  definitions 


g2n*(x-l)  =  g2n(x) 


(B2) 


and 


k2n(x-l) 


2gn(2nl ) 
(4n!) 


M2n(x) 


(B3) 


there  is  obtained 

8gn(x)  =  ~7,(2~^  f  Manfx-I)  fa'  (I)  d|  (B4) 

(4nl)  ^o 

The  changes  introduced  by  equations  (B2)  and  (B3)  are  convenient  for 
future  work.  By  equation  (B3)  the  function  M^x)  is  a  simple  fraction 
of  h3i(x-l),  the  axial  strength  distribution  necessary  to  produce  a  unit 
step  in  f2n(x).  The  functions  M2n(x)  have  been  studied  and  evaluated 
numerically  in  reference  5>  and  the  numerical  values  are  tabulated  in 
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table  I.  A  plot  of  these  functions  is  shown  in  figure  5.  These  func¬ 
tions  ,  together  with  equation  (B4),  give  a  simple  means  of  obtaining  the 
axial  strength  distribution  by  means  of  numerical  or  graphical  integra¬ 
tion  for  any  desired  f2n(x).  The  possibility  of  finding  universal 
functions  MPri  (x )  for  doing  this  is  a  consequence  of  the  fact  that  the 
body  is  cylindrical.  In  the  Karman -Moore  method  no  such  simple  rela¬ 
tionship  is  possible  because  the  body  is  of  arbitrary  shape. 
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APPENDIX  C 

ALTERNATE  METHOD  OF  OBTAINING  BODY  PRESSURE  COEFFICIENT 
OF  FOURIER  COMPONENT  COMBINATION 


A  method  of  obtaining  the  body  pressure  coefficient  involving  only 
one  numerical  integration  is  now  given.  Consider  a  disturbance  at  the 
origin  formed  by  a  bump  on  the  surface  of  the  body  as  shown  in  figure  26. 


f*  (x) 

The  nondimens ional  amplitude  of  the  bump  — — -  is  to  be  zero  except 

just  at  the  origin  where  the  area  under  the  f2n(x)  curve  is  unity,  that 
is,  the  disturbance  is  a  "delta  function".  The  amplitude  is  to  vary 
as  cos  2n 6.  Now  if  the  pressure  coefficient  in  the  wing-body  juncture 
is  -2  W2n(x)  3ue  to  the  unit  disturbance  at  the  origin,  then  the  con¬ 
tribution  to  the  pressure  coefficient  at  a  downstream  distance  x  due 
to  a  distribution  of  disturbances  from  the  origin  to  the  point  in  ques¬ 
tion  is 


^P2n  ~  ~ 


~  cos  2n Q 
V 


f2n(Ow2n(x-5)  d£ 


(Cl) 


This  integral  is  the  effect  of  all  disturbances  forward  of  the  point  and 
is  additive  to  the  effect  of  the  disturbance  exactly  at  the  point.  The 
local  value  is  that  given  by  the  Ackeret  theory,  namely, 


P  = 


(C2) 


where  is  the  local  angle  of  attack 


a7  = 


f2n(x) 


cos  2nd 


(C3) 


The  pressure  coefficient  due  to  all  disturbance  is  thus 


P2n  =  |  cos  2n6  Wx)  “  |  cos  2ne  f  W^W^x-S)  d|  (C4) 

The  advantage  of  this  method  of  obtaining  the  pressure  coefficient  at  the 
body  surface  is  that  the  pressure  coefficient  is  obtained  directly  as  a 
function  of  the  boundary  condition,  f2n(x),  without  first  determining  an 
axial  strength  function. 

The  W^Cx)  functions  introduced  to  determine  the  body  pressure  dis¬ 
tributions  are  universal  functions  that  have  been  studied  in  reference  5* 
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They  are  presented  in  table  II,  and  curves  of  them  are  presented  in  fig¬ 
ure  7*  The  values  of  WQ(x)  have  been  tabulated  by  G.  N.  Ward  in  refer¬ 
ence  11. 

Several  interesting  facts  are  illustrated  by  figure  7*  The  ampli¬ 
tude  of  the  pressure  distribution  due  to  a  unit  disturbance  of  the  type 
shown  in  figure  26  has  exactly  2n  points  of  zero  pressure,  and  these 
points  are  nearly  evenly  spaced.  The  large  number  of  points  of  zero 
pressure  for  the  higher  harmonics  makes  their  lift  contributions  much 
less  significant  than  their  pressure  contributions.  Another  interesting 
feature  of  figure  7  is  that  the  W^x)  functions  damp  much  more  rapidly 
in  the  downstream  direction  for  the  higher  harmonics  than  the  lower  har¬ 
monics.  In  fact,  for  n  =  0  the  damping  to  4  percent  of  the  maximum 
value  occurs  about  7  body  radii  downstream,  whereas  for  n  =  5  the  same 
damping  occurs  in  about  2  body  radii. 
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TABLE  I.-  VALUES  OF  Ha(x) 


X 

M0(x) 

M2(x) 

M4(x) 

M6(x) 

0 

-00 

-00 

-00 

-00 

.05 

— 

— 

-.531 

-.716 

.1 

-1.321 

-.787 

.405 

1.334 

.15 

— 

— 

— 

1.155 

•  2 

-.868 

-.166 

.902 

.735 

.25 

— 

— 

— 

.279 

.3 

— 

.131 

.831 

-.119 

•35 

— 

— 

-- M— - 

-.417 

.4o 

-.538 

.304 

.560 

-.595 

.45 

— 

— 

— 

-.675 

.50 

— 

— 

.245 

-.655 

•6 

-.378 

.450 

-.o4o 

-.417 

•7 

— 

— 

-.253 

-.070 

•75 

— 

— 

— 

.095 

.8 

-.292 

•457 

-.386 

•235 

.85 

— 

— 

— 

.340 

•9 

— 

-.434 

•  4o6 

•95 

— 

— 

— 

.433 

1.0 

-.224 

•392 

-.412 

.420 

1.1 

— 

— 

-.337 

.306 

1.2 

-.180 

.294 

-.231 

.123 

1.3 

— 

— 

-.115 

-.062 

1.4 

1 

. 

H 

5 

.188 

.000 

-.197 

1.5 

— 

— 

.095 

-.256 

1.6 

-.122 

.088 

.163 

-239 

1.7 

— 

— 

.204 

-.162 

1.8 

-.102 

.005 

.211 

-.062 

1.9 

— 

— 

.197 

.036 

2.0 

-.086 

CO 

ir\ 

0 

• 

1 

.164 

.105 

2.1 

— 

— — 

— 

.134 

2.2 

-.073 

-.100 

.069 

.121 

2.3 

— 

— 

.020 

.085 

2.4 

-.062 

-.123 

-.022 

.036 

2.5 

— 

— 

-.053 

-.013 

2.6 

-.053 

-.126 

-.073 

-.053 

2.7 

— 

— 

— 

-.063 

2.8 

-.045 

-.119 

-.080 

-.059 

2.9 

— 

— 

-.071 

-.042 

3.0 

-.038 

-.103 

-.060 

-.01? 

3.1 

— 

— 

-.038 

.005 

3-2 

-.031 

-.083 

-.019 

.022 

3-3 

— — 

— 

-.002 

.030 

3-4 

-.026 

-.061 

.011 

.028 

3-5 

— 

— 

.021 

.019 

3*6 

-.022 

-.o4o 

.027 

.008 

3.7 

— 

—  — 

.027 

-.002 

3.8 

-.018 

-.022 

.025 

-.011 

3-9 

— 

— 

.022 

-.014 

4.0 

-.016 

-.007 

.015 

-.103 

NOTE  that  Man(x)->-  7  J%  as  x—>0. 
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TABLE  II.-  VALUES  OF  W2n(x) 


X 

w0(x) 

W2(x) 

W4(x) 

W6(x) 

w8(x) 

W10(x) 

0 

0.500 

0.500 

0.500 

0.500 

0.500 

0.500 

.05 

— 

— 

— 

— 

— 

2.787 

.10 

- 

- - 

— 

2.022 

3.136 

4.459 

•15 

---- 

- - 

- — - 

- - 

3-884 

5.264 

.20 

.4319 

•7541 

1.644 

2.888 

4.176 

5.164 

•25 

— 

— 

— 

3-037 

4.021 

4.213 

•30 

— 

— 

— 

3.008 

3.486 

2.857 

•  35 

— 

— 

— 

— 

1.209 

.40 

•  3755 

.878 

1.929 

2.462 

1-702 

-.345 

.50 

— 

— 

— 

1.545 

-.215 

-2.296 

.60 

.3284 

.892 

1.428 

.508 

-1.458 

-1.986 

•70 

— 

— 

— 

-1.697 

-.461 

.80 

.2887 

.826 

.706 

-.940 

-1.116 

.980 

.85 

— 

— 

- - 

— 

1.317 

•90 

— 

— 

— 

— 

-.163 

1.345 

•95 

— 

— 

— 

— 

— 

1.107 

1.0 

.2550 

.706 

-.024 

-.974 

.634 

.676 

1.1 

— 

— 

— 

— 

•930 

-.279 

1.2 

.2262 

•  557 

-.464 

-.150 

.741 

-.803 

1-3 

— 

— 

— 

— 

.638 

1.4 

.2014 

.402 

-.558 

.457 

-.213 

-.045 

1.5 

— 

- — - 

- - 

- - 

-.479 

.399 

1.6 

1-7 

1.8 

.1801 

.255 

-.4i8 

.455 

-.514 

•475 

.251 

-.129 

-.297 

-.214 

-.005 

.154 

.167 

.061 

-.058 

-.108 

.1616 

.128 

-.142 

.073 

.022 

1*9 

2.0 

0  n 

.1454 

.026 

.081 

-.222 

.276 

d  *_L 

2.2 

O  “3 

•1313 

-.048 

.195 

-.217 

.043 

dmj 

2.4 

2-5 

2.6 

.1189 

-.097 

.196 

-•035 

-.173 

.1080 

-.122 

.120 

.106 

-•053 

2.8 

.0983 

-.128 

.025 

.098 

•  075 

.007 

3-0 

.0898 

-.122 

-.049 

-.001 

— 

.058 

3.2 

.0822 

-.104 

-.078 

-.043 

-.049 

-.0  66 

3-4 

.0754 

-.082 

-.068 

-.050 

-.030 

-.017 

3.6 

.0693 

-.058 

-•035 

-.008 

.oil 

.022 

3.8 

.0638 

-.035 

-.003 

.o46 

-- -- 

— 

4.o 

.0590 

-.017 

— 

.024 

— 
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fa)  (b)  (c) 

Figure  2.  -  Decomposition  of  wing -body  combination  into  component 

configurations. 


W,(x) 


Figure  7.— Graphical  presentation  of  W2n(x)  functions. 
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— Mach  fines 


Figure  Q.  —  Flow  field  produced  by  wing  alone . 


Figure  9.  -  Variation  of  normal  velocity  induced  at  body  surface 
by  wing  alone,  unit  body  radius . 
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Figure  //.-  Singularities  encountered  in  determining  axial  strength 

functions. 
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Figure  /5.  -  infer fei 
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Figure  16.-  Pressure  distribution  at  wing -body  juncture  of  combination  of  body  and  rectangular 

wing  using  four  and  six  Fourier  components. 


NACA  TN  2677 


4.8  erz _ 


’^1  2.4 

«3V 


| 

I  /•* 


Body 


'  Equation 

(3!  y,n*0 

\  \ 


nrQ\  \ 


Wing 


Equation 

(32)\n*i 


n-2 


n~  / 


0  1.0  2.0  3.0  4.0  5. 

y/a 

Figure  20.  -  Span  loading  of  various  Fourier  components 
acting  on  combination  of  body  and  rectangular  wing 

having  effective  chord -radius  ratio  of  4. 
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Figure  21. -Complete  span  loading  for  combination  of  body  and 
rectangular  wing  having  effective  chord-radius  ratio  of  4. 
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Figure  23.  -  Center  of  pressure  for  wing  or  control  surface  of  combination  of  body  and  rectangular 

wing. 
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Figure  25-  Approximation  of  velocity  amplitude  function 

by  series  of  step  functions. 
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6- 


Figure  2.6.  -  Fourier  component  interference  combination  with 
"delta  -  function  protruberance  at  x=0. 
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